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We calculate form factors relevant for rare B decays using moving-NRQCD for the b quark and the 
AsqTad action for the light quarks. Moving NRQCD allows us to work directly with the physical 
b quark mass and go to higher recoil momentum compared to standard NRQCD. Here, we show 
first results for the matrix elements and the operator matching coefficients. Some difficulties and 
possible ways of improvement are discussed. 
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1. Introduction 

Decays of B mesons via the flavour-changing-neutral-current transition b — > s are particularly 
sensitive to possible new-physics contributions and provide tests of the CKM mechanism at the 
loop level. Measurements of exclusive modes like B — > K*y have reached a good accuracy, and 
call for precise theoretical predictions. These are more difficult than for tree-level decays such 
as B — > TtlV, since a large set of effective electroweak operators contributes and long-distance or 
spectator effects can be important. Nevertheless, the computation of hadronic form factors in lattice 
QCD is highly desirable and complements continuum approaches. 
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Table 1: Form factors for semileptonic and radiative B decays. 



We are currently working on the calculation of the form factors listed in Table [j]. The combi- 
nation of NRQCD and improved staggered actions for heavy-light mesons has already proven very 
successful in the calculation of form factors [{j]]. In order to extend the kinematic range to high 
recoil (lower q ), we now use a moving-NRQCD (mNRQCD) action for the heavy quark. A brief 
discussion of our strategy can be found in and a new detailed account of mNRQCD will be 
given in [|p. Here, we report on the progress in the computation of matrix elements and operator 
matching coefficients achieved so far. 

2. Lattice methods 

The matrix element (F(p')\J\B(p)}, where F denotes the final pseudoscalar (P) or vector (V) 
meson and J is the relevant current in the effective electroweak operator (see Table [j]), can be 
extracted from the combination of the Euclidean 3-point correlator 

Wk^k^o^ZoHlW^ (2.1) 

y z 

with the two-point functions 

C BB (k (p) ,x ,yo) = E(*BW4»iCy))« -ik(p) ' (x ~ y) . ( 2 - 2 ) 

c FF (p',x ,yo) = £(M*)400)^ p '- (x ~ y) - (2-3) 

Here, <J>g and <J>^ are suitable interpolating fields for the initial and final meson, and /( lat > is a lattice 
version of the current, obtained by operator matching (see section |3|). For the light quarks, we 
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convert from 1 -component staggered to 4-component naive fields [QJ. Due to the use of mNRQCD 
for the b quark, the physical momenta p and q are related to the lattice momenta k( p ) and k( q ) by 

p = k {p) +Z p ym b \, 



k( q) +Z p Ym b \ 



(2.4) 



where Z p w 1 is the renormalisation of the external momentum, y = 1 / v 1 — v 2 and v is the boost 



velocity. One has p' = p — q = k( p ) 
shifted, 



The physical energy po = Eb of the B meson is also 



E B (p) =£'v(k (p) )+A v 



(2.5) 



where £ v (k( p ) ) is the unphysical energy obtained from the fit to the correlator and A v is the velocity- 
dependent energy shift. Writing t = \xq — yo\ and T = \xq — zo\, the correlators are fitted by 
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or equivalent parametrisations. Every other exponential comes with an oscillating pre-factor, as 
required by the use of naive quarks [Q. The correlator Cff receives an extra factor of 16 due 
to the trace over a 16 x 16 taste matrix, while the heavy-light correlators Cbb and Cfjb receive 
contributions from only one taste [fflj. The ground-state fit parameters are related to the matrix 
elements as follows: 
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(2.9) 
(2.10) 
(2.11) 



F = P. 

The amplitudes \[Zb, \[Zp and y/Zy in ( j2.9[ ) depend on the form of the interpolating fields <I>b, Op 
and <I>v and can be extracted from (|2.10|) and ([2.1 1[). 



3. Operator matching 

The continuum currents J must be replaced by lattice currents /( lat ) containing suitable match- 
ing coefficients to correct for the different ultraviolet behaviour of QCD and lattice mNRQCD. 
As only the high-energetic modes with E > mi, differ in the theories and GC s (mb) -C 1, matching 
coefficients can be computed perturbatively. 
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We use tadpole-improved 1-loop lattice perturbation theory. The Feynman rules are generated 
automatically []|] and diagrams are evaluated using the Monte Carlo integrator VEGAS. 

The first step is the computation of a set of heavy-quark renormalisation parameters from the 
self-energy diagrams: the zero-point energy Eq, the wavefunction renormalisation Z v , the renor- 
malisation of the mass Z m and the renormalisation of the boost velocity Z v ||]. Results for the 
full improved ff(A^ CD /ml) lattice mNRQCD action will be presented in 

Once these parameters are known, one can proceed with the calculation of matching coeffi- 
cients. For the (axial-)vector currents, these have been computed in the static limit (i.e. neglecting 
^(Agco/'Wfc) corrections in /( lat )), and the calculation including the ^(Aqcd/w^) corrections is 
underway [^]. In the following, we focus on the tensor current, which, in the continuum, is given 
by 1 
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(3.1) 



\6ti 2 ,J ' 1 ' " "~ " 2' 
We work in the static limit. At this order in the heavy-quark expansion there are two operators with 
different Dirac structure in lattice mNRQCD. For the /i = components one has 
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where m v denotes the mNRQCD field with the antiquark components set to zero. On the lat- 
tice, these operators mix under renormalisation; the one-particle irreducible vertex correction that 
contributes in the static limit is shown in Fig. [T]. 




J. 
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Figure 1: Vertex correction diagram. 



Figure 2: Matching coefficients for the tensor current. The 
subtraction point is jl = 1/a. 



Writing Jj ± = Jj\ ±7" * 2 , we obtain the lattice operator 



T (M)oe 



(\ + a s c™)J°i\ + (\ + a s cy)Jf 2 



.01' \ 



(3.3) 



The matching coefficients c± must be adjusted such that Jj 3 ^ 01 has the same one-loop matrix 
elements as the tensor operator in the continuum theory. They depend on the lattice spacing and 
contain a logarithmic ultraviolet divergence as the tensor operator is not conserved. 

'For the matching calculation, we treat the light quark as massless; in this limit the matching coefficients for the 
last two operators in Tab. |j| are equal due to chiral symmetry. 
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Results using an improved ^(Agco/^fc) mNRQCD action at am = 2.8, n = 2, the AsqTad 
action for the light quark and the Liischer-Weisz gluon action are shown in Fig. |[ The matching 
coefficient of the operator J® 1 '_, which only arises at 1-loop level, is strongly suppressed and the 
dependence on the frame velocity is found to be small for all matching coefficients. 

4. Details of the numerical calculations and first results 



In our first computations of the 3-point functions (gjj) we used the local interpolating fields 
®b(z) = q'(z)Y5b(z) and &f(x) = q'(x)r F q(x) with F P = y 5 , F v = 71,2,3. As in standard NRQCD, 
the heavy-quark Green function Gb{y,z) can be obtained by solving an initial value problem. Let 
us consider the case xq > yo > zo- Schematically, as initial value at zo we use the propagator of 
the light valence quark, 75 G q >(z,x), and then evolve the heavy-quark Green function up to the time 
slice yo where we perform the contraction with the various gamma matrices and the other light- 
quark propagator G q (x,y) = y^G] 1 {y,x)y$. This method only requires light-quark propagators with 
a fixed origin x, and, since the current is inserted only in the final contraction, allows the efficient 
simultaneous computation of arbitrary currents. 

These initial calculations were done on 400 MILC gauge configurations of size 20 3 x 64 with 
2+1 flavours of light quarks, at j8 = 6.76 and a -1 « 1.6 GeV. The light sea quark masses were 
am u = am c j = 0.007, am s = 0.05 and the light valence quark masses am u = am^ = 0.007, am s = 
0.04 (we used the AsqTad action). On each configuration, we took four different origins x, and 
additionally averaged with the time-reversed process. 

Even though we have implemented the full ^(Ag CD /mj) mNRQCD action, we only used an 
0{hQCD/mb) mNRQCD action here to save computer time. This is sufficiently accurate since we 
only considered currents in the static limit here. The heavy-quark mass was set to amj = 2.8 and 
the stability parameter was n = 2. All lattice momenta and the boost velocity were always pointing 
in 1 -direction. In this case, 21 combinations of operators/indices and final-state polarisations give 
non-zero contributions, and all the form factors listed in Table [j] can be extracted from them. 

We performed Bayesian multi-exponential fits in the two variables T and t. Gaussian priors 
for the ground state energies were taken from fit results of the corresponding two-point functions, 



with widths equal to the error from the fit result. The mNRQCD energy shift A v (see eq. (Z.5[)) was 
determined non-perturbatively from heavy-heavy meson dispersion relations (for those, the full 
mNRQCD action accurate to v\ r in heavy-heavy power counting was used). A bootstrap analysis 
was used to determine the form factors and their statistical errors. 

To give some examples, plots of the 3-point correlators ( <£> K sy b <J>jj ) and ( <£> K * sa^b <J>jj ) 
at the largest q 2 (with v = 0) are shown in Fig. || and |J The results for the tensor current in 
combination with the vector meson final state are much noisier. As expected, the statistical errors 
are seen to grow further when the recoil momentum is increased. In Fig. ||] and || we show the 
corresponding correlators at v = (0.4,0,0) and k (p) = ^(1,0,0), k( q ) = ^(2,0,0). Note that for 
the fits shown here, 4. .6 timeslices from the source/sink were skipped, so that K = 2, L = 4 (for 
the vector final state) or K = 1, L = 3 (for the pseudoscalar final state) was sufficient in ( |2.6| ). The 
results can probably be improved by extending the fitting range and using more exponentials. 

Finally, in Fig. [7] and || we show some first results for the form factors fo, f + , fx and T\, Ti. 
Note that the momentum of the meson in the final state (K or K*) was exclusively set to the very 
small values p' = or p' = ^(—1,0,0). This is made possible by the use of moving NRQCD. 
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Figure 3: Three-point correlator ( <I>a" sy^b <£>g ) at 
k,„^ =0, k (q) =0, v = 0. The fitting range is T = 14... 18 
6...(r-5). 
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Figure4: Three-point correlator (<S>k* sG\j,b Ojj )at 
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Figure 5: Three-point correlator ( <I>a" sy^b <t>g ) at 
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Figure 6: Three-point correlator (<S>k* sonb )at 
k ( P ) = f(h0,0), k (q) = f (2,0,0), v = (0.4,0,0). No 
reasonable fit was achieved yet. 
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Figure 7: The form factors /+, /o, fj (points 
for /+ and fj are offset horizontally for legibil- 
ity). The points at lowest q 2 have v = (0.4,0,0), 



f (1,0,0), k M = f (2,0,0) 
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Figure 8: The form factors T\, T%. The points at 
lowest q 2 have v = (0.2,0,0), k (p) = ^(1,0,0), 
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(2,0,0). 
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5. Conclusions 

Using moving-NRQCD and AsqTad actions, we have calculated matching coefficients for the 
heavy-light axial-, vector- and tensor currents in the static limit using 1-loop lattice perturbation 
theory, and performed first computations of form factors for rare B decays. While moving NRQCD 
significantly reduces discretisation errors at low q 2 , our initial results suffer from large statistical 
errors, overshadowing the advantages of the method. However, statistical errors do not constitute a 
fundamental obstacle and can be reduced further. The first step will be to extend the fitting range 
and include more exponentials. Then, we plan to work with random-wall sources, which were 
shown to provide considerable improvement for semileptonic decays at high recoil momentum [j|]. 
We will also use smeared interpolating fields to reduce contributions from excited states, thereby 
improving the fits. Furthermore, note that our initial computations were done with lattice momenta 
pointing in the 1 -direction only. Off-axial momenta and boost velocities will also allow lower 
values for q 2 , for example by using the final meson momentum p' = ^(— 1,— 1, 0). 

Once statistical errors are under control, we will study the dependence on the lattice spac- 
ing and the light quark masses. We also plan to include 6 ? (Aqcd / ^b) operators in the matching 
calculations. 
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